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This  work  deals  with  the  problem  of  direction  finding  using  the  matrix-pencil  approach. 
Consider  a  linear  array  of  a  sensor  and  assume  there  are  d  narrowband  sources.  Having 
collected  the  data  at  the  a  sensors,  the  problem  Is  to  estimate  the  direction  of  arrival  of 
these  sources.  The  matrix-pencil  approach  is  a  non-search  procedure,  thus  very  easy  to  use. 

In  this  study,  a  generalization  of  the  method  to  a  linear  array  of  m  identical  sensors  with 
some  arbitrary  beam  pattern  was  performed.  It  Is  also  shown  that  the  method  still  works  when 
using  different  windows.  The  only  restriction  is  that  at  least  d  elements  of  the  window  Be 
ooo-zero  to  ensure  the  validity  of  the  algorithm .  A  perturbation  analysis  due  tc  equal  sensor 
spacing  was  also  performed.  Tbs  concept  of  the  chordal  metric  was  Introduced.  It  Is  shown 
that  the  Bound  derived  on  the  chordal  metric  Is  equivalent  to  the  chordal  metric  itself.  The 
problem  of  estimating  Both  the  angular  frequencies  and  the  angles  of  arrival  of  the  sources 
was  than  posed.  It  la  proven  that  the  method  still  works;  l.e.,  the  angular  frequencies  (over 
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PREFACE 


This  work  deals  vith  the  problem  of  direction  finding  using  the 
matrix-pencil  approach.  Consider  a  linear  array  of  m  sensors  and  assume 
there  are  d  narrowband  sources,  the  signal  received  at  the  ith  sensor  is 
modeled  as 


d 

yi<t,0)  ■  £  ®l[(®|[)sij(t)  +  n^(t)  |  i»l,2, .  .  .  ,m. 


Having  collected  the  data  at  the  m  sensors,  the  problem  is  to  estimate 
the  directions  of  arrival  of  these  sources.  The  matrix-pencil  approach 
is  a  non-search  procedure,  thus  very  easy  to  use. 

In  this  study,  a  generalization  of  the  method  to  a  linear  array 
of  m  identical  sensors  vith  some  arbitrary  beam  pattern  was  performed. 

It  is  also  shown  that  the  method  still  works  when  using  different 
windows.  The  only  restriction  is  that  at  least  d  elements  of  the  window 
be  non-zero  to  ensure  the  validity  of  the  algorithm.  A  perturbation 
analysis  due  to  unequal  sensor  spacing  was  also  performed.  The  concept 
of  the  chordal  metric  was  introduced.  It  is  shown  that  the  bound  derived 
on  the  chordal  metric  is  eqivalent  to  the  chordal  metric  itself.  The 
problem  of  estimating  both  the  angular  frequencies  and  the  angles  of  ar¬ 
rival  of  the  sources  was  then  posed.  It  is  proven  that  the  method  still 
works;  i.e,  the  angular  frequencies  and  the  locations  of  the  sources 


were  estimated  using  two  matix  pencils.  The  rank  reducing  values  of 
these  matrices  is  shovn  to  contain  both  the  angles  of  arrival  and  the 
angular  frequencies  of  the  sources.  A  computer  simulation  was  performed 
each  time  to  ensure  the  effectiveness  of  the  method. 
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CHAPTER  1 


INTRODUCTION 


1.1  GOAL  OF  THE  RESEARCH 

High  resolution  direction  of  arrival  (DOA)  estimation  is  very  impor¬ 
tant  in  many  sensor  systems  such  as  radar, sonar, etc. . .  Over  the  years 
several  methods  have  been  proposed  to  solve  this  kind  of  problem.  Our 
work  is  closely  related  to  the  work  done  by  H.  OUIBRAHIM  [1].  This  ap¬ 
proach  called  the  MATRIX  PENCIL  APPROACH  addresses  the  problem  of  using 
a  passive  array  of  sensors  to  find  the  direction  of  sources  assumed  to 
be  in  the  far  field  .  The  array  is  called  passive  because  the  sources 
generate  the  signals  received  at  the  sensors.  The  received  signal  at  the 
ith  sensor  is  modeled  as 

d 

yi(t,9)  *  Z  s^Oa^)  +  nj{ t)  >1*1,2,.  .  .m  (1—1) 

k-1 

where  we  assume  the  existence  of  d  sources  and  an  array  of  m  sensors, 
ai(0k)  *s  the  relative  response  of  the  i*h  sensor  to  the  k**1 
source, 

Sfc(t)  is  the  complex  envelope  of  the  k**1  signal, 

nj(t)  is  the  additive  noise,  considered  as  the  sum  of  the  external 

and  internal  noise. 


1 


Throughout  this  work  ve  deal  with  narrowband  signals.A  problem  is 
referred  to  as  narrowband  if  the  bandwidth  of  the  impinging  signals  from 
the  sources  is  much  less  than  the  reciprocal  of  the  propagation  time  of 
the  wavefronts  across  the  array.  Hence,  given  m  measurements  collected 
at  the  sensors  we  would  like  to  estimate  the  angles  of  arrival  of  the 
sources . 

1.2  LITERATURE  SURVEY 

The  problem  of  estimating  the  location  of  sources  is  of  great  im¬ 
portance  and  has  been  approached  in  many  ways  [1-3].  Recently  several 
authors  have  suggested  a  subspace  signal  approach  [4-8].  This  approach 
is  based  on  an  eigenvalue-eigenvector  decomposition  of  the  spatial  cor¬ 
relation  matrix. This  makes  use  of  the  fact  that  there  is  a  relationship 
between  the  eigenvectors  of  the  spatial  correlation  matrix  and  the 
source  angles  of  arrival.  Moreover,  C.R.  RAO  [9]  showed  that  one  need 
only  know  the  first  few  eigenvectors  of  the  correlation  matrix.  We  now 
present  the  background  of  the  eigenstructure  approach. 

Assume  there  are  d  sources  emitting  signals  sj^t)  ;k»l,2,.  .  .  , d , 
which  are  impinging  on  a  linear  array  composed  of  m  sensors.lt  is  as¬ 
sumed  that  d  <  m. 

The  received  signal  vector  X  can  be  written  as 

X  -  A  S  ♦  N  (1-2) 
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where 

•  •  •  »xuj}“  (““I)  received  vector  signals, 

§T  -  {si,S2»  .  .  .  .sj}-  (dxl)  impinging  signals, 

*  {n^,n2,  .  .  .  ,nm}»  (mxl)  vector  noise, 

A  m  (aj  22  .  .  .  aj  }«  (mxd)  direction  matrix, 

-  (®xl)  i**1  direction  column  vector  of  A. 
ni  is  the  additive  noise  assumed  to  have  zero  mean  and  an  unknown 
variance  a2 

In  all  the  subspace  approaches  that  have  been  proposed  the  noises 
ni  are  assumed  to  be  independent  from  sensor  to  sensor  and  their  cor¬ 
relation  matrix  is  the  diagonal  matrix  a2!  where  I  is  the  identity 
matrix.  Let  the  subscript  H  denote  the  Hermitian  Transpose.  The  spatial 
covariance  matrix  is 

R  -  E[X  XH]  .  El(AS+N)(AS+N)H] 

-  BIAS  SHAH]  +  E[N  NHJ  (1-3) 

-  AE[SShJAh  ♦  a2! 

Let  S«E[S  SHJ.  Then  R  can  be  written  as 

R  *  ASA®  a2!  (1-A) 
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where 


R  is  an  (mxm)  matrix. 

Let  {Xj  <  X2  <  X3  <,  .  .  .  , <  Xm  }  be  the  set  of  eigenvalues  of  R. 

Let  {Yi,V2,V3,  .  .  .  , Vj,  }  be  the  set  of  the  corresponding  eigenvectors. 
If  S  is  non  singular  and  with  the  assumption  that  m  >  d,ve  can  show  that 

1)  the  minimum  eigenvalue  of  R  is  with  multiplicity  (m-d),i.e, 

. V^nin-*2- 

2)  the  eigenvectors  associated  with  the  minimum  eigenvalue, 

Yd+2'  ¥d+3»  *  •  •  »Ym»  are  orthogonal  to  the  space  spanned  by  the 
columns  of  A.  This  can  be  written  as 

{Yd+l»Yd+2'Yd+3’  ■  ’  *  »Ym^  JL  fal»a2,a3'  *  *  *  »ad^  (1-5) 

where 

_L  ,  denotes  orthogonality, 
a^  ■  i1*1  column  of  A. 

This  algorithm  can  be  summarized  as  follows: 

1)  determine  the  number  of  sources  d  from  the  multiplicity  of  Xmjn. 

2)  the  orthogonality  relation  (1-5)  between  the  direction  vectors 
of  the  impinging  sources  and  the  eigenvectors  corresponding  to  Xm^n 
yields  the  directions  of  arrival  of  the  sources.  Ve  just  have  to 
"search"  for  those  directions  vectors  that  are  orthogonal  to  the  eigen¬ 
vectors  corresponding  to  Xm^n.  For  this  reason  these  methods  are  called 
search  procedures.  They  assume  the  eigenvalues  and  the  eigenvectors  to 
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be  perfectly  known.  However, in  practice  this  is  not  always  true.  Ve  then 
need  to  perform  some  kind  of  optimal  estimation.  Such  a  procedure  was 
developed  by  LIGGET,LAVLEY  and  BARTLETT  [10-12].  In  such  procedures  some 
hypothesis  testing  is  introduced.  One  difficulty  is  the  subjective 
judgement  required  to  set  the  thresholds.  An  approach,  which  is  consid¬ 
ered  as  one  of  the  best,  is  the  Aikake  Information  Criterion  (AIC), since 
it  does  not  require  any  objective  judgement  on  the  thresholds.  Another 
approach  is  the  Minimum  Description  Length  (MDL)  approach.  But  one  of 
the  most  promising  techniques  is  Multiple  Signal  Classification  (MUSIC) 
proposed  by  SCHMIDT  [5].  This  algorithm  provides  asymptotically  unbiased 
estimates  of 

1)  number  of  signals, 

2)  directions  of  arrival, 

3)  strengths  and  crosscorrelation  among  the  directional  waveforms, 

4)  polarizations 

5)  strength  of  noise/interference. 

More  recently,  other  methods  are  being  developed.  Some  of  them  are  non¬ 
search  procedures.  These  approaches  have  very  importants  advantages  over 
search  procedures.  A.  PAULRAJ.R.  ROY  and  T.  KAILATH  [8],  in  their  ap¬ 
proach  known  as  ESPRIT,  have  shown  that  their  algorithm 

1)  does  not  require  knowledge  of  the  array  geometry  and  element 

characteristics  (directional  pattern, gain/phase) , 
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2)  is  computationally  much  less  complex  because  it  does  not  use  the 


search  procedure, 

3)  does  not  require  a  calibration  of  the  array, therefore  eliminat¬ 
ing  the  need  for  the  associated  storage  of  the  array  manifold  which  can 
be  very  large  for  multidimensional  problems, 

4)  simultaneously  estimates  the  number  of  sources  and  DOA's  . 


1.3  OUTLINE  OP  THE  WORK 

H.  OUIBRAHIM  [1]  proposed  a  generalization  of  the  ESPRIT  method. 
This  method  consists  of  applying  an  operator  to  the  received  signals  in 
order  to  form  a  matrix  pencil  M-AN.  The  rank  reducing  values  of  A  are 
shown  to  contain  the  information  needed  to  estimate  the  DOA's. 

The  pencil  theorem,  presented  in  chapter  2,  establishes  the  rela¬ 
tionship  between  the  rank  reducing  values  of  A  and  the  functional  form 
f(4j)  generated  by  the  operator  applied  to  the  measurements. 

In  chapter  3  a  generalization  of  the  method  to  arbitrary  but 
identical  beam  patterns  is  presented.  Both  the  cases  of  deterministic 
signals  and  zero-mean  random  signals  are  considered. 

Previously,  only  rectangular  windows  have  been  applied.  In  chapter 
4  it  is  shown  that  the  method  still  works  using  different  windows.  As  in 
chapter  3,  the  cases  of  deterministic  signals  and  zero-mean  random  sig¬ 
nals  are  considered.  A  comparison  of  the  different  windows  is  obtained 
by  means  of  a  computer  simulation. 
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In  chapter  S  a  perturbation  analysis  for  the  case  of  deterministic 
signals  is  performed.  The  concept  of  CHORDAL  METRIC  introduced  by 
STEVART  [13]  is  used.  The  chordal  metric  is  a  very  good  measure  of  the 
perturbation  between  the  perturbed  eigenvalue  and  the  true  one.  A  bound 
is  derived  which  is  shown  to  be  effictive  by  means  of  a  computer  simula¬ 
tion. 

Chapter  6  is  devoted  to  a  new  technique  for  the  simultaneous 
estimation  of  the  angular  frequencies  and  the  angle  of  arrival  of  d 
sources  assumed  to  be  in  the  far  field.  The  technique  makes  use  of  the 
decomposition  of  two  (2)  matrix  pencils. 

Finally,  a  summary  and  some  suggestions  for  future  work  are  given 


in  chapter  7. 


CHAPTER  2 


REVIEW  OF  MATRIX  PENCIL 


Our  problem  is  the  estimation  of  the  angles  of  arrival  of  d  sources 
given  measurements  collected  at  the  m  sensors.  The  expression  for  the 
received  signal  presented  in  equation  (1-1)  of  chapter  1  shows  that 
the  measurements  are  linear  combination  of  d  exponentials  whose  ex¬ 
ponents  j+k  ;  k»l,2,...,d,  contain  the  information  needed  to  determine 
the  locations  of  the  sources.  Specifically, 

♦it  -  «/c  D  sinO^)  ;k«l,2, ....  ,d.  (2-1) 

The  pencil  theorem  establishes  the  relationship  between  the  rank 
reducing  values  of  X  and  the  functional  form  f(Aj)  generated  by  the 
measurements. 

2.1  PENCIL  THEOREM 


Denote  by  C  the  field  of  all  complex  numbers. Consider  two  matrices 
M  and  N  of  size  (kxp).The  set 

(  M-AN  ;  X  e  C  J  is  said  to  be  a  pencil. 

The  matrices  M  and  N  are  required  to  have  the  following  decomposition s 
M  -  E  P 
N  -  E  D  F 

where 
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E  is  a  (kxd)  matrix  and  k  >  d 
P  is  a  (dxp)  matrix  and  p  >  d 
0  is  a  (dxd)  diagonal  matrix. 

Theorem 

* 

If  M  and  N  are  two  matrices  which  have  the  decompositions  cited 
above  and  if  E,F  and  are  all  of  rank  d,  then  the  rank  of  the  matrix 
pencil  M-XN  is  decreased  by  1  whenever 

xi  *  <dii>-1  f  i-l»2,...,d.  (2-2) 


Proof 

Since  N-EF  and  N-EDF, 

M-XN  -  EF-XEDF 
-  E(I-XD)F. 

Thus 

rank  (M-XN)  «  rank(E(I-XD)F) 

■  min{rank(E) , rank(F) ,rank(I-XD)} . 
However,  by  assumption 

rank(E)«rank(F)«d 

and 

rank(I-XO)  is  of  rank  d  as  long  as 

l“xidii«0* 

If  1-Xjd^i»0  which  implies  that  xi*(dii)“*  ,rank(I-XD)»d-l. 
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Therefore,  the  rank(I-XD)  is  reduced  by  1  whenever 

Xi«(dii)  *  ;  i-l>2, . . . . ,d. 

In  our  work, the  matrix  0  has  all  its  entries  of  the  form 

eJ^i  . 

Thus,  the  rank  reducing  values  of  Xj  are 

Xi  -  (eJ+i)-1  -  e-3+i  ;i-l,2, _ ,d. 

2.2  EVALUATION  OP  THE  RANK  REDUCING  VALUES 

We  have  assumed  previously  the  existence  of  d  sources  and  a  linear 
array  composed  of  m  sensors  vith  the  condition  that  d  <  m/2.  We  then 
formed  two  matrices  M  and  N  of  size  (m-d)xd.  We  see  that  tvo  cases  may 
occur.  If  m-2d,  M  and  N  are  tvo  square  matrices.  The  set  of  the 
generalized  eigenvalues  of  the  pencil  M-XN  is  defined  to  be  the  set  of 
all  elements  Xj  such  that 

det(M-XiN)-0. 

When  the  generalized  eigenvalues  are  distinct, the  rank  of  M-XN  is 
reduced  by  1  whenever  X  equals  one  of  these  values.  In  the  case  where 
d  <  m/2,  M  and  N  are  non  square  matrices.Det(M-XjN)  no  longer  exists 
since  the  pencil  is  not  square.  For  this  reason  we  have  to"make"  the 
pencil  matrix  a  square  one.  This  can  be  done  by  premultiplying  the 
pencil  M-XN  by  either  M®  or  N®. 


10 


V«  obtain 


MH(M-XN)  -mhm-xmhn 


or 

Nh(M-XN)  «NaM-XNaN. 


Ma(M-XN)  and  Na(M-XN)  are  both  square  matrices  (dxd).  Notice  that 
Ma<  M- XN ) - ( EF ) a<  EF- XEDF ) « FaEaEF- XFHEHEDF 
«FaEaE(I-XD)F 

and 

Na ( M- XN )«( EDF ) a <  EF- XEDF ) «  FHDHEHEF-  XF^E^DF 
«FaDaEaE(I-XD)F. 

Both  equations  have  the  decompositions  required  by  the  pencil  theorem 
since 

FaEaE  and  F  are  of  rank  d 


and 

FaDHEHE  and  F  are  of  rank  d. 

Because  (I-XD)  arises  in  all  of  these  equations ,  ve  can  say  that  the 
generalized  eigenvalues  of  MH(M-XN)  and  Na(M-XN)  are  identical  to 
those  obtained  for  the  case  vhere  M  and  N  were  square  matrices. 


il 


CHAPTER  3 


GENERALIZATION  TO  ARBITRARY  BUT  IDENTICAL  BEAM  PATTERN 

The  moving  vindov  developed  by  H.  Ouibrahim  {1]  was  shown  to  be 
a  non  search  procedure.  It  was  applied  for  azimuth  only  DOA  (direction 
of  arrival)  estimation  of  far  field  point  sources.  A  generalization  of 
this  approch  to  arbitrary  but  identical  beam  pattern  is  presented  here. 

3.1  DETERMINISTIC  CASE 

Assume  we  have  a  linear  array  composed  of  tn  identical  sensors 
vith  uniform  spacing  D.  Assume  there  are  d<  m/2  narrovband  sources  lo¬ 
cated  at  azimuthal  angles  0jt,k-lf2, . . .  ,d,  which  are  impinging  on  the  ar¬ 
ray  as  planar  wavefronts  and  emitting  signals  whose  complex  envelopes 
are  denoted  by  sjjt)  ,  k-1,2  ,3,...,d.  The  received  signal  at  the  ith 
sensor  is  modeled  as 

d 

t )aj(0j^)+nj ( t )  j i*l * 2 , . . . . m,  (3—1) 

k*l 


where 

ai(®k)  is  the  relative  response  of  the  i ^  sensor  to  the  kt^1 
source, 

nj(t)  is  the  additive  noise  assumed  to  be  zero-mean  Gaussian. 
ai(9k)  can  be  written  as 


«! 


« 
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ai<ek)-a<©k)  eJ<i-1><w/c)  D  sin<0  >  ;i-l,2,...,m,  (3-2) 

vhere 

w  s center  frequency  of  each  of  the  spatial  sources , 
c  : speed  of  propagation  of  the  plane  waves, 
a(6)  :  beam  pattern  of  each  sensor. 

If  ve  let  ♦k  «  «  D  sin(0|t)/c,  yj(t,0)  can  be  rewritten  as 

d 

yiCt,©)-!  sk(t)a(9k)e^i-1>*k  +  njCt)  , i-1,2, . . . ,m.  (3-3) 

Taking  the  expected  value  of  equation  (3-3),  we  get 

d 

x1(t,0)-E[yi(t,0)]  -  I  sk(t)a(9k)eJ<i-1)^k  ; i-1,2, . . . ,m.  (3-4) 

k-1 

Define  the  rectangular  window 

(  1  i  i$n$N 

RN(n)"< 

(  0  ;elsewhere 

»  Xj ( t , 0)  , 

component 
2n 


Given  the  number  of  sources  d  and  the  m  averaged  data  points 
ve  create  (d+1)  vectors  Xn  ;  n-1,2,.  .  .  ,d+l  , where  the  i1*1 
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xn+i-l  ( 1  *  ®)Rtn-d(  * )  * 


n»l f  2  i  *  * • j d+1 
1$  i  <  m-d. 


Specifically,  if  the  argument  (t,0)  is  omitted  for  simplicity, 

XXT  ■  (X1 »x2 *••••* xm-d) 

X?*^  *  (x2»x3'  *  •  *  *  ,xm-d+l^ 


XdT  ■  (xd»xd+l» ' ‘ * ,xm-l^ 

Sl*lT  *  (xd+l»xd+2* • *  *  * »xm^ • 

The  matrix  pencil  H-XN  is  then  formed  vhere 


t  t 

I  I 

2l  *2 

!  I 

1  l 


t  t 


*2  b 


l  1 


Letting  a^aaCd^)  and  omitting  the  argument  t  in  s^Ct),  Xj  can  be  writ¬ 


ten  as 
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xi 

x2 

alsl  +il2s2  +  •  1/ . +  adsd 

aiSjeJ*!  +  a2S2eJ"  + . +  *dsdeJ  ° 

5i  - 

• 

• 

- 

• 

• 

xa-d 

alsl®^  +  . +adsde3  (n-d--l)#d 

M«J*1  12.2*!*2 


•  ad 

•  ade 


J*d 


^(n-d-D^jdn-d-l)^ 


^(B-d-l)#d 


In  general,  X„  is  given  by 


xn 

aisi®^ ♦  •  • 

.  .  +  adsdeJ<a*l^d 

*n  - 

xn+l 

• 

• 

» 

m 

aisje^n^l  +  .  .  .  . 

• 

• 

• 

.  .  +  adsdeJn*d 

• 

xm-d+n-l 

a1s1e^n”^+n-^^  ♦ 

.  .  .+  adsde^n-d+n-^^d 

a}eJn+l 


a2«j(n-l)42 

a2e^n^2 


.  .  .  a,jej 
.  .  .  ade^d 


a1eJ<a”d+n-2)*l  a2e^m_d+n_2)*2  .  .  .  adeJ(n-d+n-2>^d 


Xjj  can  also  be  written  as: 


Xjj  -  A  *(n-l)  BS 


(3-5) 


and  the  matrices  M  and  N  become 


t  t  f 


ABS  AB*S  AB^S . AB*(d*1)s 


i  i  i 


N. 


t  T  t 


AB*S  AB*2s  ABpS . Afl^S 


■i.  4.  -t 


Factoring  out  AB  in  M  and  AB*  in  N  ,  we  get 


M  .  AB 


t  t  t 


S  *S 


♦<d-l>< 


ill 


•J 
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Let  F  be  the  matrix 

■  t  I1  t  j 

F»  S  .  *<d-l)s 

(  i  i  4> 

and  E  the  matrix 

E»AB . 

Ve  see  that  M  and  N  have  the  decompositions 

N  -  E  F 
N  -  E  *  F. 

The  matrix  pencil  then  becomes 

M-AN  «  EF-XE*F 

=.  E(I-A*)F  (3-6) 

which  satisfies  the  requirements  of  the  pencil  theorem.  Since  E*AB, 
rank(E)«min(rank(A) , rank(B) ) .  But  A  is  of  rank  d  as  long  as  the  direc- 
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tions  of  arrival  of  the  signals  are  distinct  and  X  <  d/2.  B  is  of  rank  d 
as  long  as  a^>a<d^)  , i»l,2, . . . ,d,  are  different  from  zero.  F  is  of  rank 
d  even  in  presence  of  coherent  sources. 

Hence,  the  values  of  X  for  which  the  rank  of  M-XN  decreases  by  1  are 
given  by 

X^  »  e“3^k  ;k»l,2, . . . ,d. 

The  angles  of  arrival  are  given  by 

0^  »  sin-l(jcln(Xfc)/cod)  ;  k=l,2,...,d. 

3.2  ZERO- MEAN  RANDOM  CASE. 

As  before  the  signal  'received  at  the  ith  sensor  is  modeled  as 
follovs; 

d 

yi(t,e)«E  sjt(t)ai( 0|^)  +n^ (t)  ;i»l,2, . . .  .m. 
k-1 

Analogous  to  the  previous  section,  we  form  (d+1)  vectors  In  in-1,2 . ■ 

where 

*  (yi ( t , 9)  ,y2 ( t ,  0) ,  •  .  ♦  • 

I2T  -  {y2(t»0)*y3<t»®)»-  •  •  •  'ym-d+l^'®)) 


Idll“  {yd+l(t’®)»yd+2<t’0>*  •  •  •  •»ym^*®^* 
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Define  the  inner  product 


«h,k  ■  <  Ih-Ik  >  ■  EI  Sc"Xh  I- 
Ve  then  form  the  matrices  and  as  follows 


Ml 


"1 


■1,1  “1,2  •  •  •  •  ®l,d 
“2,1  “2,2  •  •  •  •  “2,d 


“d,l 

“d,2  • 

•  *  '  “d,d  j 

“2,1 

“2,2 

....  ®2 ,  d 

“3,1 

“3,2 

....  m3fd 

•  •  • 


“d+1,1  “d+1,2 


•  “d+ld 


The  vector  In  can  be  decomposed  as 


In- 


yn 

yn+i 

• 

ais^e^n-^^^l  + 
alsle^n^-  +  •  •  • 

•  •  •  +  adsde^*_1)*d 

.  .  .  +  adSde-J"^ 

+ 

• 

• 

ym-d+n-l 

alsl*^B1_^+n~^^^  + 

nn+l 


nn+m-d-l 


Let  Njj  denote  the  noise  vector.  ^  can  be  written  as 
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L«t  S  be  the  vector 

ST  ■  {s^ > S2 »••••» S(j  ) 
Then  Yjj  can  be  revritten  as 


In  •  A  ♦<n-1>  BS  +  Njj  (3-8) 

vhere 

1  1  .  •  •  1 

eJ+1  eJ*2  ...  e^+d 

•  «  • 

•  •  • 

A  - 

gj (m-d-l)^!  gj(m-d-l)^2  ,  .  .  ;j(n-d-l)*d 
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Ve  knov  chat 


®h,k  -  CXh’L^  *  <  AB#^h"1^S+Nh  ,  ABf<k-1)s+llk  > 

-E[  ^.^1  -  EI(AB4^k-1)s+Nk)H(AB*<h-1)s+Nh)J.  (3-9) 

Assuaing  the  signals  and  noise  to  be  statistically  independent,  ve  can 
vrite 

®h,k  -ElSH*H(k-1>BHAHAB*<h-1>S]  +  E [ N®jcNh ] .  (3-10) 

Suppose  the  noise  components  are  statistically  independent  Gaussian  ran¬ 
dom  variables  with  zero-mean  and  variance  <?-.  Then 


E 


0 

(m-d)o^ 


;k#h 

jk-h 


(3-11) 


and  mhfk  is 


mh,k* 


E|SH*H(k-1)BHAHAB#(h-l>S]  ;  kjth 

E[SH*H^k-1)BHAHAB*<h-1)s]  +  (m-d)o2  ;  k-h. 


A  more  convenient  expression  for  is  nov  derived.  Note  that 


1  e--j*l  . 
1  e-j+2  . 


1  e-J*d  . 


^(m-d-l)^ 

e-j(m-d-l)A2 


;_j(m_d-l)*d 


1 

e^l 


1 

e-J*2 


1 

e-1  ♦d 


ej (m-d-1) ^  gj (m-d-l)$2 


ej (m-d-1) 
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Multiplying  #B(k~l)A®A  by  ve  get 


*H(k-l)AHA4(h-l), 


<m-d)eJ<h-^h  F,1e-j(^-l)*leJ<h-1)*2 

Flje-JOc-l'-  — -  21  —  — * 


■l)^(h-l)*l  <5id)eJ<h-k>*2 


•Fdle 

•pd2« 


-j(lc-l>*leJ(h-l>*d 

-j(k-l)*2ej(h-l)*d 


Flde:J<k-l)+deJ(h-l)h  F2d«-J<k:l>W<h-l>*2.  .  .  (m-d)eKh-k)#d 


aiSi(«-d)eJ  +  a2a2^21e”^k“^^l  e^h-l)^  + 

. ♦  adsdFdle-Klc-l)h  eiO'-D+d 

eJ  +  a2S2(m-d)eJ^"^^^2  + 

. •*■  adsdFd2e"^*t_^^^2  eJ  (^-1)^ 


alsl/l2e”^k”*^  eJ(h-l)^i  +  a2S2F2de_^^_^^d  e3(k-l)$2 
. +  adsd(m-d)ej^~*c^d 


and  becomes 

a*s*[  ajSj(a-d)eJ (h-lt)^i  +  a2s2^21e~^k“^^l  ej (h-1) ^  + 

. .  adsd^dla”^  (k-l)+i  eJ  (h-l)^d  ] 

+ 

ajs^l  als1^12®”^ +  a2S2(m-d)e3(k-lt)+2  + 
•  •  •  •  +  adsd^d2®~^k“^^2  ] 


adsdl  als1^12®-^ (^-1) Aj  ej(h-l)+^  +  a2S2F2de~^ (k-l)$d  ej(k-l)$2 
. +  a^SjjCm-dJeJ^^^d  J. 


Noting  that  ej(h-k.)$j  can  ^e  written  as 

•J(h-k)4i  .  e-j<fc-l>*i  eJ<h-D*i  for  all  i-l,2,.  .  . 
and 

Pii  -  (m-d)  for  all  i-l,2,...,d, 


ve  obtain 


d  d 


sHBH|B(k-l)AHA*(h-l)BS  «  E  E  Ppqsjspajape-^-1^  e^h_1^p. 


q-1  p-1 


If  ve  let 
S 


Ets*sn  ] 


Pq  “iaqaP 


»pq 


aqV 


then 


d  d 


E(SH*B<k"1)BHAHAB*<h-1)s]-  E  E  SpqapqFpqe-j  eJ(h_1)*p. 

q-1  p-1 


mh>k^ecomes 


“h 


•A 


E[sH^<lc-l)BH4HAB,(h-l)s1  ;h* 

E|sH*H(k-l)BHAHAB#<h-l)si  +  (m-d)o2  ;h-k 


(3-13) 


r  E  E  SpqapqPpqe-^^-1)^  ej<h-l)*p 
q-1  p-1 


;h*k 


(3-14) 


d  d 

E  E  SpqapqFpqe-^"1^  ei(h~1>*p  +  (m-d)o2  ;h-k. 

q-1  p-1 


Let  I  be  the  identity  matrix  and  1^  the  matrix  defined  as  follovs 
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Il 


0  10  0.  . 

0  0  10.  . 

0  0  0  1.  . 


0  0  0  0.  . 

0  0  0  0.  . 


Finally  ve,  define  the  matrices  M  and  N  as  follovs 
H  -  Mi  -  (m-d)ff2  I 

N  -  Nx  -  (m-d)ff2  l1.  (3-15) 

The  matrix  pencil  is 

M-XN  .  (M1-(m-d)ff2I)-X(N1-(m-d)«»2I1). 

Define  the  matrices  U,  V  and  4  as  follovs  : 


1  1  ...  1 

eJ^l  eJ^2  .  .  .  eJ^d 

•  •  • 

*  •  • 

U  ■  .  . 

^(d-l)*!  ej(d-l)+2  .  .  .  ej<d-l)*d 


sllallFll  s12a12F12  •  •  •  •  SldaldFld 

s21a21p21  s22a22F22  •  *  •  •  s2da2dp2d 

V  - 

sdladlFdl  sd2ad2Fd2  •  •  •  •  sddaddFdd 
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■  ■ 

*3*2 

0 

♦  * 

0 

It  can  be  shown  that  H  and  N  have  the  following  decompositions  : 

M  -  U  V  UH 

N  -  U  V  4s  UH.  (3-16) 

Hence,  as  required  by  the  pencil  theorem,  the  matrix  decomposition  of 
the  pencil  M-XN  is  given  by 

M-XN  .  (UVUH)-X(UV*BUa) 

«  UV(I-X#)UH.  (3-17) 

The  matrices  UV  and  UH  are  of  rank  d  as  long  as  the  directions  of  ar¬ 
rival  of  the  signals  are  distinct.  Therefore  the  values  of  X  for  which 
the  rank  of  the  pencil  M-XN  ■  ^(I-X*0)!]®  is  decreased  by  1  are  given  by 

X4  -  eJ*i  ;i-l,2,...,d.  (3-18) 

The  angles  of  arrival  are  given  by 

0i  -  sin-1(-jln(Xi)/«D)  ji-1,2,.  .  .  ,d.  (3-19) 


) 
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CHAPTER  4 
WINDOWS 


In  chapter  3  the  rectangular  vindov  was  the  logical  and  obvious 
choice.  It  wa a  used  to  form  the  sequence  Xn+i-l(l» 5)®m-d(*)*  In  this 
chapter  ve  show  that  any  shaped  vindov  vill  vork;  i.e,  the  angles  of  ar¬ 
rival  are  obtained  from  the  matrix  decomposition  of  a  matrix  pencil.  A 
computer  simulation  and  a  comparative  performance  of  the  different 
vindovs  used  are  also  presented. 

4.1  DETERMINISTIC  CASE  : 


Consider  a  linear  array  composed  of  m  identical  sensors  vith 
uniform  spacing  D.  Assume  there  are  d  <.  m/2  sources  emitting  signals 
vhose  complex  envelopes  are  denoted  by  s^Ct)  ;  k»l,2, .  .  .  ,d. 

As  before  the  received  signal  at  the  i'*1  sensor  is  modeled  as 

d 

y^t.O)  -  Z  sk(t)  a^e^)  +  n^( t)  ;  i«l,2,...,m  (4-1) 

k-1 


vhere 

a^Q^)  is  the  relative  response  of  the  i**1  sensor  to  the  kth 
source, 

nj(t)  is  the  additive  noise  assumed  to  be  zero-mean  Gaussian, 
a^ (6)  can  be  vritten  as 


1 
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where 


a^ej-aOJe^1-1^ 


M«D/ c)  sin(6). 

For  simplicity,  the  arguments  (t,Q)  in  y^(t,8)  and  t  in  both  S|<(t)  and 
nj(t)  are  dropped. 

Taking  the  expected  value  of  (1)  we  get 
d 

xi  ■  E[y±]  >  I  S|(  aCQfcJej^-^^k  ; i-1,2, . . . ,m.  (*-2) 

k-1 

Consider  the  sequence  xn+i-lvm-d(*)  ;n«l,2, . . . , d+1  and  i«l,2, . . . . ,m-d, 
where  is  the  value  of  the  window  of  width  (m-d)  evaluated  at  the 

point  (i).  Ve  then  form  (d+1)  vectors  Xjj  where 

xnwm-d(l) 

xn+lwm-d(2) 

• 

xn+m-d-l^m-d ( ) 

If,  for  simplicity,  we  let  cj  •VB_<|(i),  £n  can  be  written  as 
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■  ' 

*  * 

C1 

xn 

c2 

xn+l 

0 

• 

• 

0 

• 

cm-d 

xn+m-d-l 

However,  it  is  shown  in  chapter  3  that 


xn+l 


-  A#<n_1>BS 


(4-3) 


xn+m-d-l 


where 


1 

e-J 


A 


^(n-d-l)^  gj(a-d-l)^2 


gj(m-d-l)*d 
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ST  ■  {sj_  S2 . Sjj}  • 

Let  C  be  the  matrix  given  by 


Then  Xjj  can  be  written  as 

Xn-CA*<n-1>BS. 
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This  expression  holds  for  various  choices  of  the  c*  corresponding  to 
different  vindovs. 

The  matrices  M  and  N  which  are  formed  by 


’  t  t 

I  I 

. .  .T  ' 
1 

'  t 

1 

t 

1 

. . .  T 

1 

H  « 

? 1  X2  •• 

;  N  - 

h 

*3 

?d+l 

1  1 
.  ±  ± 

1 

l  . 

1 

.  i 

1 

1 

1 

become 


III  I 

CABS  CABtS  CAB^S . CAB*(d“1>S 


I  I 


N 


III  I 

CABtS  CABt^S  CAB^S . CAB^S 


If  ve  factor  out  CAB  from  H  and  CABt  from  N,  we  get 
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N  -  C  A  B  ♦  F, 
the  matrix  pencil  becomes 


(4-5) 


M-XN  -  CABF-XCAB*F 

-  CAB(I-X#)F.  (4-6) 

Let  B  denote  the  matrix  CAB.  Then 

M-XN  -  E(I-X#)F  (4-7) 

which  satisfies  the  requirements  of  the  pencil  theorem.  Note  that 
rank(E)*min{rank(A) ,rank(B),rank(C)} . 

We  have  seen  that  A  is  of  rank  d  as  long  as  the  directions  of  arrival 
are  distinct  and  the  separation  0  is  less  than  X/2.  B  is  of  rank  0  as 
seen  earlier.  When  choosing  the  elements  of  the  matrix  C,  it  is  neces¬ 
sary  that  at  least  d  of  the  diagonal  elements  be  non  zero.  This  vill 
ensure  that  rank(C)  >  d.  Therefore  the  rank  of  the  pencil  M-XN  is 
decreased  by  1  whenever 

Xt  -  eJ+i  ; i-1,2, .  .  .  ,d.  (4-8) 

4.2  ZERO-MEAN  RANDOM  CASE  ; 

Again  the  signal  received  at  the  i**1  sensor  is  modeled  as 
d 

/i(t»9)  *  E  ^k(  t)  &i(0jr)  +  n^(t)  t  i-“l > 2 , . . . , m.  (4-9) 

k-1 

yi(t,8)  -  xi(t,0)  +  nj(t)  ;i=l,2, . . . ,m.  (4-10) 
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Consider  the  sequence 


yn+i-l^’^m-d^)  ;n«l,2f . . .  ,d+l 

and  1<  i  <  m-d. 


For  simplicity,  the  arguments  (t,9)  and  t  are  dropped,  (d+1)  vectors  Yjj 
are  then  formed  vhere 


Xn 


ynvm-d^^ 

yn+l^m-d^^ 


yn+m-d-lVd(m-d> 


As  before,  let  Cj^W^^Cj).  Then 


P 

» 

r 

C1 

xn 

C1 

nn 

c2 

xn+l 

c2 

nn+l 

0 

• 

0 

• 

• 

• 

+ 

• 

• 

• 

• 

• 

• 

0 

. 

0 

• 

cm-d 

xn+m-d-l 

cm-d 

nn+m-d-l 

. 

L  J 

Let  C  be  the  matrix 


i 
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Noting  that 


An 

xn+l 


xn+m-d-l 


A*(n-1>BS, 


In  can  be  written  as 

Yn  *  CAB*^n_1^S  +  CNn  (4-11) 

where 

Nn  =  {nn  nn+i . nn+m-d-l  }  • 

This  expression  holds  for  various  choices  of  c^  corresponding  to  dif 
ferent  windows. 

Define  the  inner  product 

mh,k  =  <Ih>Ik>  -  MlkHIhl- 

Define  the  matrices  and  as  follows 


®1»1  ®1»2 . ml’d 

®2  *  1  ®2 » 2  *  *  *  •  *  ®2  ’  d 


«1  « 


®d»  1  ®d*2 . ®d*d 


4 


®2»1  ®2»2  ®2*d 

®3»1  ®3»2  ®3»d 

»  •  • 

N*  - 

•  •  • 

•  •  • 

®d+l»l'md+l*2 . ®d+l»d 


Since  Yjj  *  CAB*(n_1^S  +  CI^,  it  follows  that 
»h,lc  ■  E[Ika.Yhl 

.  E[(Sa*H<k-1>BflAHCH  +  NkHCH)(CAB*<h-1>S  +  C^)]. 

Assuming  that  the  signals  and  noise  are  statistically  independent, 
®h,k  -  EtSH*H<k-1>BHAHCHCAB*(h-1)s  ]  +E[NkHCHCNh) ] . 

Let  the  noise  components  be  statistically  independent  Gaussian  random 
variables  with  zero-mean  and  variance  <P- .  Then 


E[NkHCHCNh)] 


'  0 

m-d 

<P-{  E  |cA  |2  ) 
i-1 


;h#k 

;h-k. 


(4-12) 


% 


I 
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Assuming  all  c j's  to  ba  real.  Icil2”ci2'  Therefore  ,  CHC  is  the  matrix 


CHC* 


c2* 


0 


cm-d 


It  follows  that 


CHCA. 


cl2 

1  1  ...  1 

c22 

eJ  ^1  e  J  ^2  .  .  .  e-i  ^d 

0 

• 

•  •  • 

•  •  • 

0 

•  •  • 

•  .  .  •  • 

cm-d2 

ej(m-d-l)#i  ej(m-d-l)*2  .  .  .  eJ<m-d*1>^d 

Cl“  , 

C2^eJ  ^1 


=2**1 +2 


.  .  .  ;m  d2ej(«-<l-l)+d 


m-d 


If  we  define  Fpq  »  £  ej(i-l)($p  -  $q)  f  then  AHCHCA  becomes 


i-1 
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F11  f21  •  •  •  *  Fdl 

f12  f22  •  •  •  •  Fd2 


AhChCA  - 


•  • 
Fld  F2d 


Following  the  same  procedure  as  in  section  3.2, 


it  can  be  shown  that 


E[sHtB(k-l)BB*BcHCAB«<h-1>Sl.  £  I  Sp,ap(,Fp,e-l O1'1  > ♦<! 

q-1  P-1 


where 


Spq-E[SqSpl 

apq-ajap 


,"idc,2  ,J<i-D<V  *q>  . 
i-1 


Therefore, 


t 


mh,k* 


E[SH^^lt“1^BHAHCHCAB*(h“1)SJ 


E(SH*B^k'1^BHAHCHCAB*(h_1^Sl  + 


m-d 

J-i,  t  |ci|2  ) 

i-1 


;h*k 

;h-k 


4.5 


‘  C  \  SnWtqt-Hl- !)♦,  ej(h-l)*p 
q-1  p-1 

1  £  spqapqFpqe~j(k-1)*q  ♦  <*2(  z  IcJ2 

lq-1  p-1  i-1 


Let  I  be  the  identity  matrix  and  1^  the  matrix  defined  as 


;h*k 

(4-13) 

;h-k. 


0  10  0 
0  0  10 
0  0  0  1 


O'  0  0  6 
0  0  0  0 


0 

0 

0 


1 

0 


Also,  define  the  matrices  M  and  N  as  follows  : 


m-d 

M-Mi  -  (  I  cj2)  a2  I 
i-1 


N-Nj  -  (  Z  Ci2)  a2  li  . 
i-1 


The  matrix  pencil  becomes 


m-d  m-d 

M-AN  -[Mi  -  (  Z  cj2)  a2  IJ-AfNi  -  (  E  cj2)  a2  I^. 
i-1  i=l 
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If  we  define  the  matrices  U,  V  and  #  by 

1  1  .  .  .  1 

ej*l  ei+2  ...  eJ*d 

•  •  ♦ 

•  •  • 

U  - 

ej(d-l)*i  ;j(d~l)^2  .  .  .  e^-^+d 

sllallFll  s12a12F12  ♦  *  *  *  SidaldFld 

s21a21F21  s22a22F22  •  •  •  •  s2da2dF2d 
#  •  • 

•  *  • 

V  - 

,  •  • 

•  *  * 

sdladlFdl  sd2ad2Fd2  •  •  •  •  sddaddFdd 
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it  can  be  shovn  that  M  and  N  have  the  folloving  decompositions 
M  -  U  V  UB 
N  *  U  V  f3  UH. 

The  matrices  UV  and  UH  are  of  rank  d  as  required  by  the  pencil  theorem. 
Therefore  the  rank  of  the  pencil  M-XN  -  UVCI-X*®)!!8  is  decreased  when¬ 
ever 

Xj  «  e^i  ;i*l,2, . . . ,d. 

The  angles  of  arrival  are  given  by 

©i  ■  sin“^(-jcln(Xi)/ftiD)  ;i-l,2,.  .  .  ,d.  (4-14) 

Hence,  theoretically  any  shaped  vindov  would  give  the  same  directions  of 
arrival  of  the  sources. 


4.3  COMPUTER  SIMULATION 

In  this  section  the  comparative  performance  of  the  rectangular, 
Hamming,  Hanning  and  Blackman  windows  is  evaluated  by  means  of  a  com¬ 
puter  simulation. 

The  different  windovs  are  defined  as  follows: 


Rectangular 


RN<n) 


1  ;  0  <  n  <  N 
0  ;  elsewhere 
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"1 


I 


Hanning 


%(n) 


.5(l-cos(2nn/N))  ;  0  <  n  <  N 
0  ;  elsewhere 


Hamming 


RN(n) 


.  54- .  46cos  ( 2  Jtn/N )  ;  0  <  n  <  N 
0  ;  elsevhere 


Blackman 


RN(n) 


.42-.5cos(2rm/N)+.08cos(4rm/N)  ;  0  <  n  <  N 
0  ;  elsevhere. 


The  scenario  used  for  this  simulation  consisted  of  two  coherent  sources 
(d«2)  which  are  incident  on  a  linear  array  consisting  of  eight  sensors 
(m-8).  The  sources  are  assumed  to  be  located  at  9j*18°  and  02*22°  . 

The  received  signal  at  the  i1^  sensor  was  modeled  as 


yj(t,0)  «  Z  sj^t)  a(0jt)ej -t-nj(t)  ;  i*l,2,...,m, 
i-1 


where 
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♦k  -  «D/c  sinCe^). 

In  this  simulation  the  sensors  were  chosen  to  be  omnidirectional.  Con¬ 
sequently, 

a(0k)  -  1 

Also  the  noise  vas  simulated  to  be  white  Gaussian  with  zero-mean  and 
unit  variance.  The  sensors  were  positioned  at  half  wavelength  apart  such 
that  coD/c  *  n  .  Finally,  the  complex  envelopes  were  selected  to  be  » 
s 2  -s  where  s  is  a  constant.  The  signal  to  noise  ratio  is  defined  as 

PS 

SNR  -  - 
Pn 


The  cases  considered  in  this  simulation  are  shown  in  table  (4-1) 


Table 

(4-1) 

SNR 

|s| 

30  dB 

22.36 

10  dB 

2.24 

In  this  simulation,  equation  (4-15)  which  gives  the  angles  of  arrival, 
was  not  used  because  it  assumes  that  magnitude  of  X  is  unity.  Because 
this  was  not  the  case  in  actual  practice  due  to  numerical  inaccuracies, 
the  method  did  not  perform  well.  To  overcome  this  situation  the  follow¬ 
ing  approach  was  used. 
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Let  X  -  a  +jb  vhere  a  and  b  are  tvo  real  numbers.  Using  an  exponential 
notation  X  can  be  written  as 

X»  {a2+b2}^exp{jtan“l(b/a)} . 

Prom  equation  (4-14)  X»exp{j4}.  Thus,  ignoring  the  magnitude  of  X 

4  »  tan“l(b/a) 

and  the  angles  of  arrival  are  given  by 

0  m  sin-1  {(ctan-l(b/a))/«D} . 


The  results  of  the  simulation  are  shown  in  tables  (4-2)  and  (4-3). 

.  Table  (4-2) 


Rectangular 

Hanning 

Hamming 

Blackman 

SNR 

1  ©1 

02  | 

01  02 

1  ©1 

02  | 

1  ©1 

©2 

30  dB  | 

|  17.982 

|  21.995  |  17.882  {  21.948  |  17.887 

{  21.999  18.217 

22.418 

10  dB 

18.050 

|  22.142 

17.295  |  23.672 

18.294 

|  23.324 

|  14.529 

|  21.065 

Mean  of  6^  and  @2 
(  500  snapshots/run  10  runs) 
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Rectangular 
®1  02 

Hanning 

®1  02 

Hamming 

01  ®2 

Blackman 
®1  e2 

■  30  dB  |  0.018  | 

0.072  |  0.126 

|  0.078  |  0.177  | 

0.090  |  0.138  |  0.317 

10  dB  |  4.385  | 

4.960  |  11.418 

|  14.654  |  7.843  | 

9.552  |  20.639  |  4.595 

CHAPTER  5 


PERTURBATION  ANALYSIS 


In  this  chapter  ve  want  to  investigate  the  behavior  of  the 
moving  vindov  in  the  presence  of  perturbations  due  to  sensor  spacing.  Ve 
restrict  ourselves  to  first  order  perturbations.  STEWART  has  shown  [13] 
that  a  good  measure  of  this  perturbation  is  the  chordal  metric  which  is 
introduced  next. 

5.1  CHORDAL  METRIC 


Let  C  denote  the  field  of  all  complex  numbers.  Consider  two 
matrices  M  and  N  and  let  X  be  an  eigenvalue  of 

M  x  «  X  N  x.  (5-1) 

x  is  called  right  the  eigenvector  of  equation  (5-1).  Also  let  jr  be  a 
left  eigenvector  of  the  matrix  pencil.  £  satisfies 

M  -  X  N.  (5-2) 

For  convenience,  x  and  £  are  usually  normalized.  Thus,  we  set 

11*11  -  1  «d  |  Ie|  |  -  1. 

Ve  also  Introduce  the  Euclidean  matrix  norm  defined  as 

I  Ml  -  sup  |  |Mx  1 1 . 

Ve  are  interested  in  the  generalized  eigenvalue  problem 
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(5-3) 


M  x  -  X  N  x. 


vhere 


M-M+AH-M  +  E 

N-N  +  AN-N  +  F 


(5-4) 


Let 

a  -  2®Mx 

(5-5) 

and 

0  -  £%x 

(5-6) 

From 

equation  (5-1)  it  follovs 

X  «  o/0. 

(5-7) 

Stevart  [13]  shoved  that  small  perturbations 

in  E  and  F  result 

o+£HEx+0( e2 )  a' +0( e2 ) 

0*2HFx+O(e2)  0'  +0(e2) 

(5-8) 

vhere 

0(c2) 

lim  -  »  0. 

e-*0  e 

Define  the  chordal  metric  as 

i  >■->.  i 

X(X,X)  - - . 

(5-9) 

\  1+ |x j2  1+|X|2 
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Geometrically  X(X,  X)  is  half  the  length  of  the  chord  connecting  X  and  X 
vhen  they  have  been  projected  in  the  usual  vay  onto  a  Rieman  sphere  of 
unit  radius.  The  maximum  value  of  the  chordal  metric  is  unity. 

With  this  definition  it  vas  shown  [13]  that 

X(X,X>  <  c/y  +  (He2)  (5-10) 

where 

e  «<!  |  |E 1 12  +  |  |F  1 12  (5-11) 

r  -J  a2  +  02  (5-12) 

and 

a  and  0  have  been  defined  earlier. 

5.2  APPLICATION  TO  THE  MOVING  WINDOW 


Again  assume  a  linear  array  composed  of  m  identical  sensors 
spaced  at  D+ADi  where  dD^*0.  Assume  there  are  d  narrowband  sources.  The 
received  signal  at  the  i**1  sensor  is  modeled  as 

d 

yi(t>0)—  1  Sjf( t )a<  ( 8i» )  +  (t)  «  i=l»2» .  • .  >ui  (5—13) 

k-1 
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where  the  """  denotes  the  response  of  the  perturbed  array, 

nj(t)  is  the  additive  noise  assumed  to  be  zero-mean  gaussian, 
and 

ai(Ok)  is  the  perturbed  relative  response  of  the  i1*1  sensor  to  the 
kth  source. 

Taking  the  expected  value  of  equation  (5-13)  ve  have 
d 

Xi(t,0)»  E  sv(t)a< (0k)  1 i— 1 , 2 ,  •  • » ,m .  (5—14) 

k-1 

Note  that 

ai(9)  -  a(0)exp{^i~1^D^<o/c^  sin(0)-t-j(ft>/c)dDisin(0)j 

-  a(0)eJ(i-1)D<“/c)  sin(0)  ej(«/c)dDisin(0). 

To  a  first  order  approximation 

ej(w/c)dDisin(0)  a  l+j (oj/c)dD1sin(0)  -  1+j (2ndDi/5)sin(0) 
where  6  is  the  va.'elength  of  the  signal  wavefront. 

Thus  aj(0)  can  be  vritten  as 

aj(0)  -  a(0)eJ<i-1>D<“/c>  sin<0>+j(2xdDi/5)  eJ ( i-l)D(«/c)  sin(0)sin(e)a(0) . 
Equation  (5-14)  becomes 
d 

Xi(t,0).  E  (a(0k)sk(t)(ei<i-1)D<«/c>  sin(0k) 
k-1 

♦jUn^/S)  e3(i-l)D(«/c)  sin(0k)sin(ek)a(ek))>  (5-16) 

For  simplicity,  denote  a^-atO^).  Then 


(5-15) 
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d 

XiCt.e).  I  aksk(t)(eJ<i"1)D(tt/c)  sinOfc) 
k-1 

d 

+j(2nADi/6)  I  aksk(t)(eJ<i'1)Dw/c  sin<ek)sin(0k). 
k-1 


(5-17) 


Notice  that  the  first  part  of  equation  (5-17)  is  just  the  non-perturbed 
quantity  xj  which  appeared  in  chapters  3  and  4.  Dropping  the  argument 
(t,9)  in  equation  (5-17)  ,  it  can  be  vritten  as 


xj-  xj  +  toq  -  xj  +  e*.  (5-18) 

(d+1)  vectors  X„  are  then  formed  vhere  Xjj  is  given  by 

«•* 

X  / 

In  *  {  xn  xn+j . xn+m-d-l  )  • 

Xfj  can  be  vritten  as 


xn 

*n+l 

• 

*  * 

en 

en+l 

• 

In  « 

• 

• 

xn+m-d-l 

+ 

• 

• 

en+m-d-l 

Xn  *  ^ 


(5-19) 


where  E^'  -  {  en  en+j . en+m-d-l  )•  In  chapter  3  and  4  it  is 

shovn  that  X„  can  be  expressed  as 


Xn-AB*<n-1)s 


(5-20) 
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vhere  B  is  nov  the  identity  matrix  and  A,  ♦  and  S  are  : 


1  1  ...  1 

eJ^l  e^^2  .  .  .  e^d 

•  •  • 

•  •  • 

A  *  •  •  • 

gj(m-d-l)$igj(m-d-l)+2  .  .  .  *j(m-d-l)*d 


sT  -  {  sj  S2 . } . 

The  matrices  M  and  N  and  then  formed  in  the  usual  vay.  Specifically, 


‘  T 

1 

T 

1 

t 

1 

■  t 

i 

t 

1 

T 

I 

N  - 

Si 

S2 

. Xd 

;N  - 

S2 

S3  •  • 

•  •  -Sd+i 

1 

.  1 

1 

1 

1 

1 

1 

.  i 

1 

4, 

1 

1 
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To  apply  equations  (5-10)  through  (5-12),  one  has  to  get  the  Euclidean 
noras  of  E  and  F. 

The  Euclidean  matrix  norm  is  defined  as 

I  lE  1 1  -  sup  |  |Ex  1 1 . 

Mil  l-i 


Hovever 

| |Ex| |-{(Ex)H(Ex)}%  -  (xhEhEx}%.  (5-23) 

For  simplicity,  let  D«EEE.  Notice  that  maximizing  ||Ex||  is  the  same  as 
maximizing  the  quadratic  form  xEDx.  The  maximization  of  xEDx  involves 
obtaining  the  eigenvector  of 

Dx  *  Xx 

vhich  corresponds  to  the  largest  eigenvalue  of  D.  Then 
max(xflDx)-XE|naxxHx.  (5-24) 

Since  xHx«l . 

max(x0Dx)-XEmax 

and 

MElM^max)1*  (5-25) 

where  XEm|)X  is  the  largest  eigenvalue  of  EHE.  Similarly, 

i  lF  I  l*(^max^  (5-26) 

vhere  XEmax  is  the  largest  eigenvalue  of  FHF. 

It  follows  that 
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c  -1  |  |E|  I2  7  ||f||2'  .  \  +  XFmax  .  (5-27) 

Ve  know  that  D«E®E  is  a  symmetric  matrix.  Let  djj  be  the  ii**1  element  on 
the  diagonal  of  D  where  djj  is  real.  Define  the  TRACE  of  D  as  the  sum  of 
all  the  elements  on  the  digonal  of  D  which  also  equals  the  sum  of  all 
eigenvalues  of  D.  Hence  XEmax  <  trace(D)-sum(dn> .  To  obtain  the  trace, 
the  diagonal  elements  of  D«EEE  must  be  computed.  Recall  that  E  is  the 
matrix 


I 

Therefore, 

I 


«1  *  e2  •  •  •  •  «d 

e2  e3  ....  ed+1 


em-d  em-d+l 


E". 


L 

_  ★ 

* 

eU 

e2* 

e2 

e3 

• 

«d* 

• 

ed+l 

*m-l 


em-d  * 
em-d+l 


*m-l 


D  is  then  the  matrix 
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el  e2 
e2  e3 


ed 

ed+l 


* 

•  •  •  •  “n-d  ^ 

•  •  •  •  em-d+l 


D  - 


em-d  em-d+l 


em-l 


ed+l 


ea-l 


* 


Since  ve  are  interested  in  the  trace  of  D,  only  its  diagonal  elements 
are  enumerated  belov,  where  use  is  made  of  the  fact  that  e}*0. 


D 


m-d 

Z 

i-2 


m-d+1 

£  l*i  I2 

i-2 


m-1 

?  I«i 

i-m-d 
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trace(EHE)  -  |e£  I2  +  |«3  I2  - +  l«m-dr 

+|«2l2  +  l«3 I2  + - +  lem-d I2  +  lem-d+ll2 

♦  le3  I2  lem-d  I2  +  lem-d+l  I2  +  lem-d+2l2 


+  M2  + 


led+ll2 


l«m-l  l: 


For  i*l,an  element  ej  of  E  can  be  written  as 


«i  *  j 


£  ajj  e-3(i_*)*k  sin(9jc) 
k-1 


Note  that 


2it 

l«il  «  -  IdDil 

4 


£  a^e^1-^^  sin(9k)  | 

k-1 


2k  d 

<  -  I^Di!  E  I  sk  ak  eJ^-^^k  sinO^) 

k-1 


2n  d 

<  -  |  ADi  |  E  laksk 


$ 


2n  d 

£  l^max I  z  lak  sk  I*  (5-30) 

S  k-1 

Hence, 

2n  d 

trace(EHE)%  £[(m-d-l)+(a-d)(d-l)l%  -  l^maxl  E  laksk  I 

&  k-1 


trace(EHE)lfe  <ld(m-d)-l]%  -  l^maxl  1  lak  sk  !•  (5-31) 


Using  exactly  the  same  procedure,  it  can  be  shovn  that 


2n  d 

trace(FHF)%  £[d(m-d)J%  -  |ADmaxl  z  lak  sk  I* 

S  k-1 


Since 


(5-32) 


xEmax  <trace(EHE) 
and  XFmax  <trace(FHF), 


C-JllEll2  ♦  l|F||2  -  J  xEnax  *  X*B 


2n  i - -  d 

<  -  |dDmax|  J  2d(m-d)-l  (  E  |aksk|  ). 

6  k-1 


Finally,  with  reference  to  equations  (5-10)  through  (5-12),  the  bound  on 
the  chordal  metric  becomes 
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2ji  I  i  d 

-  |ADmax|  <1  2d(m-d)-l)  (  E  |aksk|  ) 

5  k-1 

X(X,X)  S  - 

•I  (£HMx)2  +  (£HNx)2 


(5-33) 


where  jrH  and  x  are,  respectively,  the  right  and  left  eigenvectors  cor¬ 
responding  to  X. 

A  second  approach  for  obtaining  a  bound  on  the  chordal  metric 
makes  use  of  the  Frobenius  norm  which  is  defined  for  an  (mxn)  matrix  A 
by 


I  I A  |  |f  -{  ”  E  |aij|2  . 
i-1  j-1 


(5-34) 


It  can  be  shown  [14]  that  the  Euclidean  norm  is  always  less  than  or 
equal  to  the  Frobenius  norm.  Specifically, 


I 


l|A||f  >  | |A | | 


Recall  that  E  is  the  matrix  defined  as 


el 

... 

•  ed 

e2 

e^  ... 

• 

• 

•  ed+l 

em-d 

em-d+l  '  ‘  * 

•  em-l 
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where  is  equal  to  zero. 

Also,  from  equation  (5-30), 

2  it  d 

M  ±  -  l^maxl  1  lak  sk  I* 

S  k-1 

It  follovs  that 

d 

|  |  E  |  |f  <  {  <2n/5)2  |  |2  <<m-d)d-l)(  E  |ak  sjtl)2}1*.  (5-35) 

k-1 

Similarly,  we  have 

d 

||F||f  <  {  (2rc/5)2  |  ADmax  |2  ((o-d)d)(  E  |ak  SjJ)2}*  (5-36) 

k-1 

Therefore, 

«  "I  llEll2  *  l|F||2'  <  l  l|E||t2.||F||£2 

2n  I - -  d 

<  I ^max  I  2d(m-d)-l)  (  E  I ak.sk  I  )• 

6  k-1 


The  bound  on  the  chordal  metric  then  becomes 


X(X,X)  < 


(5-37) 


—  I^maxl  J  2d(m-d)-l)  <  E  |aksk|  ) 

Notice  that  this  bound  is  exactly  the  same  as  the  one  obtained  while 
using  the  Euclidean  norm. 


5.3  COMPUTER  SIMULATION 

To  evaluate  the  usefulness  of  the  derived  bound  for  the  chordal 
metric  a  computer  simulation  was  performed.  The  model  used  consisted  of 
two  deterministic  sources  (d«2)  with  constant  and  equal  complex  en¬ 
velopes  incident  on  a  linear  array  consisting  of  eight  sensors  (m«8). 

The  sources  are  assumed  to  be  located  at  and  02-22°.  For 

simplicity,  the  case  of  omnidirectional  sensors  vas  assumed.  In  the 
simulation  the  case  of  perfect  sensor  spacing  vas  first  considered.  500 
snapshsts  were  used  to  obtain  the  matrices  M  and  N.  The  process  vas 
repeated  10  times  and  the  results  averaged  to  obtain  nominal  values  for 
Xj,  jo  and  jrj  ;i-l,2.  A  pertubation  AD  -.001D  vas  then  introduced  and 
the  procedure  used  in  the  unperturbed  case  was  repeated.  The  unperturbed 
signal  received  at  the  i^  sensor  was  modeled  as 

d 

yt(  t,  0)-  Z  s^tJeJ*1"1)"70  Dsin^)  +  ni(t)  ;i=i,2, . . . ,m 
k-1 


whereas  the  model  for  the  perturbed  signal  was 


d 

yWt,®)-  £  slc(t)e-i<<i-1)D+aDi>w/c  sin<®k>  +  ni(t)  ;i.i,2, . . .  ,m. 

~  k-1 

nj(t)  vas  simulated  as  white  Gaussian  noise  with  zero-mean  and  variance 
o^«l.  D  was  assumed  to  be  equal  to  S/2  (half  the  wavelength). 

Because  s^«S2*s  ,  the  signal  to  noise  ratio  is  is  defined  to  be 
SNR-2  |s|2. 

The  cases  considered  are  listed  in  table  (5-1) 

Table  5-1 

1  SNR  j  ]s |  | 

30  dB  |  22.36  | 

|  10  dB  |  2.24  | 

The  computed  results  are  shown  in  table  (5-2) 
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Observe  that  the  bounds  for  the  chordal  metric  are  of  the  same  order  of 
magnitude  as  the  chordal  metric  itself. 


CHAPTER  6 


SIMULTANEOUS  ESTIMATION  OP  DIRECTIONS  OP  ARRIVAL 
AND  ANGULAR  FREQUENCIES 

In  the  previous  chapters  the  sources  vere  assumed  to  be  emitting 
at  a  common  angular  frequency  u  with  complex  envelopes  denoted  by  s^t) 
;k*l,2,.  .  .  ,d.  In  this  section  ve  show  that  the  matrix  pencil  method 
still  works  for  the  case  where  the  signals  have  different  center  fre¬ 
quencies.  The  angles  of  arrival  and  the  angular  freqencies  are 
estimated  simultaneously. 

6.1  DETERMINISTIC  CASE 

Consider  a  linear  array  of  m  identical  sensors  uniformly  spaced 
at  a  distance  D.  Assume  there  are  d  <  m/2  narrowband  sources  located  at 
azimuthal  angles  6^  ;k-l,2(.  .  .  ,d,  which  are  impinging  on  the  array  as 
planar  wavefronts  and  emitting  signals  whose  complex  envelopes  are 
denoted  by  s^Ct)  ;k«l,2,.  .  . ,d.  The  first  sensor  is  followed  by  an 
equally  spaced  tapped  delay  line  consisting  of  m  taps  with  successive 
delays  of  T  seconds  (see  fig  6-1). 

With  reference  to  the  first  sensor,  the  signal  received  at  the  ith 
sensor  is  modeled  as 

d 

y^( t, 0)  -  r  a(0k)sk(t)  eJ(i-1>D<wk/c>sin^9k>  +  n^t)  ji-1,2 . m.  (6-1) 

k-1 
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Th«  signal  received  at  the  (h+l)c**  delay  line  tap  is 
d 

yi(t-hT,0)  -  E  a(6k)sit(t)eJhT“k  +  ni(t-hT)  ;h-0,l . m-1.  (6-2) 

“  k-1 

Assume  the  noise  to  be  a  zero-mean  stationary  Gaussian  process.  Then 
taking  the  expected  values  o£  equation  (6-1),  ve  get 

d 

x1(t,e)-Elyi(tf0)l-  E  a(©it)sk(t)  eJ<1-1>D<“k/c>sin<ek>;i-l,2,.  .  .,m.  (6-3) 

k-1 


Similarly,  the  expected  value  equation  (6-2)  is 


zh( t-hT, 0)-E[yi( t-hT, 9) ] 


d 

*  E  a 
k-1 


(0k)sk(t)eJhT“k  ;h-0, 1, . 


. ,  m— 1 . 


(6-4) 


Let  ♦jt«(ujfD/c)sin(0jt);  k-1, 7, .  .  .,d.  Equation  (6-3)  becomes 
d 

xi(t,0)«  E  a(0jl)S)c(t)  eJ^'^^k  ;i»l,2,.  .  .,m.  (6-5) 

”  k-1 

For  simplicity,  the  arguments  (t,0)  in  Xi(t,9),  (t-hT, 9)  in  Zh(t-hT,9) 
and  (t)  in  sj^t)  are  dropped.  Denote  a)c-a(9ic).  Equations  (6-5)  and  (6-4) 

become 


x^-  E  a^Sjj  eJ^"^)^k  ;  i— 1 , 2 , . 

k-1 


.  ,m, 


(6-6) 
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and 


% 


■  £  ajfSuc 

k-1 


4**“k  jh-0,1,.  .  . ,m-l. 


(6-7) 


We  then  form  (d+1)  vectors  and  (d+1)  vectors  Zjj,  vhere 

?n  *  (xn  xn+l  •  •  •  •  xn+m-d-l^»  n»l,2,.  .  .  ,d+l, 
and 

ZrT  -  (Zr-1  zr  *  *  *  *  zr+m-d-2l»  r»l,2,.  .  . ,d+l. 
hi  can  be  expressed  as 


xn 

xn+l 

• 

• 

M 

aiSjeJ^11”1^  +  .  . 

a^sieJn^l  +  .... 

.  +  a^sje^11"^-^^ 

.  +  adsde-in*d 

• 

• 

xm-d+n-l 

alSlej<“-d+n-2)*l  +. 

.  +  adsdeJ<m-d+n-2>*d 

hi 


a-i  (n_^) 

aieJnh 


a2ej<n-l)+2 

a2eJn^2 


ade^  (n-l)+d 
ade^ 


a1eJ<a-d+n-2>h  a2eJ(m"d+n'-2)^?  .  .  .  adej<m-d+n-2>*d 


S1 

s2 


sd 


Xjj  can  also  be  written  as: 
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(6-8) 


Xn  -  A  *<n-l)  b§. 

Following  the  same  procedure,  it  can  be  shown  that  In  can  be  written  as 
Zn  -  AjY^"1)  BS  (6-9) 

where  B  and  S  have  been  defined  previously  and  A^  and  Y  are  the  matrices 

1  ...  1 

eJT“2  .  .  .  eJTfifc 

•  •  ? 
gj (m-d-l)T<i>i  gj (m-d-l)Ta>2.  .  .  (m-d-DTuy 


gjToh 

eJTw2 


’ejTwd 


Four  (4)  matrices  M,  N,  P  and  Q  are  then  formed,  where 


r  t  t 

1  1 

. .  t  ■ 
1 

■  t 

1 

t 

1 

t 

•  •  •  1 

1 

M  - 

h.  h.  •• 

••  xd 

;  N  » 

-2 

*3 

xd+i 

1  1 

.  a  a 

1 

4. 

1 

.  A 

1 

A 

1 

A 
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t  t  t 


N  -  AB«  S  *S 


*(d-l)S 


4-4-4. 


Let  F  be  the  matrix 


rtt 


F«  S 


*(d-l)S 


4-4-4- 


and  E  the  matrix 


E-AB. 


H  and  N  have  the  decompositions 


H.EF 


N  .  E  ♦  F. 


The  matrix  pencil  M-XN  can  then  be  written  as 

M-XN  -  EF-XE#F 


(6-10) 


E(I-X*)F 


(6-11) 


which  satisfies  the  requirements  of  the  pencil  theorem.  Hence,  the 
values  of  X  for  which  the  rank  of  M-XN  decreases  by  1  are  given  by 
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^  •  6  ^ { k««  1|2|  •  •  •  f 


(6-12) 


d. 

The  same  procedure  is  now  applied  for  the  decomposition  of  the  matrices 


P  and  Q.  It  follows 

‘  t  t  t  t 

P  -  AjBS  AiBTS  AjB^S . A1B^d~1>S 

.  ill  i 

t  T  T  t 

Q=  ARBI'S  AjB^S  A1BT3S . A^BI^S 

i  4  i  ^ 


Factoring  out  AjB  in  P  and  AjBY  in  Q,  we  get 


'  T  j  T  T 

P  «A]B  S  VS  f2S .  f(d-l)s 

i  i  A 
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t  t  t 


Q  -AjBY  S  YS  Y2S 


Y  <  d — 1  >  s 


4-  4-  i 


Define  the  matrix  G  as 


t  t  t 


G-  S  YS  Y2S .  Y<d_1>S 


14-  4- 


Let  E^  be  the  matrix  AjB.  P  and  Q  have  the  following  expressions, 


P  -  Ej_  G 
Q  «  Ej  Y  G. 

The  matrix  pencil  P-5Q  becomes 

P-8Q-  E^-SEiYS 
-  E1(I-SY)G, 


(6-13) 


(6-14) 


which  satisfies  the  requirements  of  the  pencil  theorem.  Hence,  the  rank 


reducing  numbers  of  the  pencil  P-5Q  are  given  by 

■  e"^^  ;k*l,2,  —  ,d. 


(6-15) 
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Note  Chat  aquation  (6-15)  gives  us  an  estimate  of  w^.  This  estimate  is 
used  in  equation  (6-12)  to  get  an  estimate  of  9^.  Therefore,  this  meth¬ 
od,  simultaneously  estimates  the  angular  frequencies  and  the  angles  of 
arrival  of  the  sources;  i.e 

»i  -  jln(5i)/T 

0i  -  sin-1((jln(Xi))/(«i  D)) ;  i-1,2 . d.  (6-16) 


6-2  ZERO-rfKAN  RANDOM  CASE 


In  this  section  ve  assume  that  the  complex  envelopes  of  the 
emitted  signals  are  stationary  random  processes  vith  zero-mean.  It  is 
shovn  that  the  matrix  pencil  method  still  vorks;i.e,  the  angular  fre¬ 
quencies  and  the  angles  of  arrival  of  the  sources  are  estimated  simulta¬ 
neously  using  matrix  decompositions  of  tvo  (2)  matrix  pencils.  Again, 
let  the  received  signal  at  the  ith  sensor  be 

d 

yi(t,0)  -  E  a(0u)sk(t)  ej(i-l)D(Vc)sin<9k>  +  nt(t)  si-1,2,.  .  .,m.  (6-17) 

k-1 

The  signal  received  at  the  (h+l)'*1  delay  line  tap  is 
d 

yx( t-hT, 0)  -  E  a(0lc)sk(t)eJhT“k  «•  n^t-hT)  ;h-0,l . m-1.  (6-18) 

k-1 
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Let  4ka(ukD/c)sin(0k);  k«l,2,. 


.  .  ,d.  Equation  (6-17)  becomes 


d 

yi<t,0)-  Z  a(0!c)sk(t)  +14(0;  i-1 ,2, .  .  .,m.  (6-19) 

k»l 


For  simplicity,  the  arguments  (t,0)  in  Xi(t,0),  (t-hT,0)  in  zh(t-hT,0) 
and  (t)  in  sk(t)  and  in  n^(t)  are  dropped,  denote  ak»a(0jt).  Equations 
(6-19)  and  (6-18)  become 


d 

Yi«  1  *ksk  eJ<1_1)*k  +  nt  ;  i»l,2,.  .  .,m,.  (6-20) 

k»l 


and 


d 

zh”  E  akske^hT<*k  +  nh  jh-0,1,.  .  . ,m-l.  (6-21) 

k-1 

Analogous  to  the  previous  section  (d+1)  vectors  Yj,  and  (d+1)  vectors  Zj 
are  formed,  where 

Xn  *  (Yn  Yn+l  •  •  «  •  Yn+m-d-1^'  n»l,2,.  .  . ,d+l, 

and 

T 

Xn  *  fzn-l  zn  •  •  '  *  zn+m-d-2^;  n"l»2,.  .  .  ,d+l. 

Define  the  inner  products  mk>k  and  win  to  be 

mh,k-<  Ih’Ik  >» 

and 
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Four  matrices  Mj,  N^,  Pj  and  Qj  are  then  formed  as  follows 


and 


w2, 1  v2,2  •  •  •  •  w2,d 

w3, 1  w3,2  •  •  •  •  w3,d 
•  •  • 

Qj  » 

•  •  • 

•  •  • 

vd+l , 1  vd+l,2  •  «  •  •  vd+lfd 

Following  the  same  procedure  as  in  chapters  3  and  4,  it  can  be  shown 
that  Yjj  and  Zr  can  be  expressed  as 

Yn=AB*(n"1)s  ♦  Njj 

and 

Zr.A1BY<r-1)s  +  Nr . 

A,  Alf  B,  ♦,  Y,  and  S  have  been  defined  earlier,  N„  and  Nr  are  given  by 


T 

Nn  ■  {  nn  nn+j  ....  nn+(n_d_j  } , 

and 

5r^  *  f  nr  nr+l  •  •  •  •  nr+m-d-l  i • 

For  simplicity,  denote  k-1,2,.  .  .,d.  The  matrix  Y  then  becomes 
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We  know  that 


®h,k  *  <Jh»2k>  *  <  AB*<h"1h+Nh  ,  AB^-^S+Nfc  > 

-El  Yjj®. YjjJ  -  El(AB*<k-1>S+Nk)a(AB*(h-1>S+Nh)l.  (6-22) 

Assuming  the  signals  and  noise  to  be  statistically  independent,  ve  can 
write 

“h,k  -E[Sfl*H<k-1>BHAHAB*<h-1)s]  f  E[NHuNh].  (6-23) 

Suppose  the  noise  components  are  statistically  independent  Gaussian  ran¬ 
dom  variables  with  zero-mean  and  variance  a2.  Then 


EISkgShl' 


0  ;k*h 

(m-d)a2  ;k-h 


(6-24) 


and  m^'it  is 


mh,k* 


EISH*H(k-l)BHAHAB#(h-l)s]  •  k*h 

(6-25) 

EjgHfH(k-l)BHAHABt(h-l)si  +  (m-d)o2  ;  k-h. 


Ve  obtain 


SHBB#B(k-l)AHA#(h-l)BS  ,  z 

q-1 


1  FpqsqspaqaPe"'i<k"1)^q  eJ(h~1)*P 

p-1 


where 


•J 
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/jdej(i-:L)(*p  -  *q)  . 
i-1 


Hence, 


ElsH*H(k-1)BHAHAB*(h-1)s]-  l  I  SpqapqFpqe-^-1)^  •Jth-D+p 

q-1  p-1 


where 


Spq**E[  SqSp  ] 

apq*aqap  • 

f  ElSH*H<k-1)BHAHAB*<h-1>S]  ;hjtk 

,k  l  E[SH*H<k-1)BeAHAB*<h-1)s]  +  (m-d)o2  ;h-k 


(6-26) 


1  Z  SPqapqFpqe  J(k_1>+q  ej(h_1>*p  ;h*k 

q-1  p-1 
d  d 

E  £  spqapqFpqe~^k_1^q  ej(h_1>*p  +  (m-d)c2  ;h-k. 

lq«l  p-1 


(6-27) 


Let  I  be  the  Identity  matrix  and  1^  the  matrix  defined  by 
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0  10  0. 

0  0  10. 

0  0  0  1. 


.  0 

.  0 

.  0 


•  •  •  • 

0  0  0  0 

0  0  0  0 


1 

0 


Finally,  we  define  the  matrices  H  and  N  as  follovs 


M  «  -  (m-d)o^  I 

N  -  Nx  -  (m-d)o2  Ij.  (6-28) 

The  matrix  pencil  is 

M-AN  -  (M1-(m-d)<j2I)-X(N1-(m-d)a2I1). 

Define  the  matrices  U,  V  and  ♦  as  follovs  : 


« 


I 


1  1 

eJ*l  eJ*2 


U 


*  •  • 

ej(d-D*l  ;j(d-l)+2  .  .  .  e^d'D+d 


sllallFll  s12a12F12  •  •  •  •  sldaldFld 
s21a21F21  s22a22F22  •  •  •  •  s2da2dF2d 


sdladlFdl  sd2ad2Fd2  •  •  •  •  sddaddFdd 


« 
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i 


It  can  be  shown  that  M  and  N  have  the  following  decompositions 


H  .  U  V  U“ 

N  •  U  V  *8  Ua. 


(6-29) 


Hence,  as  required  by  the  pencil  theorem,  the  matrix  decomposition  of 
the  pencil  M-XN  is  given  by 


M-XN  >  (UVUH)-X(UV*HUH) 
-  UVCI-X^JU0. 


(6-30) 


Therefore,  the  values  of  X  for  which  the  rank  of  the  pencil  M-XN  »  UV(I- 
X*H)Ufl  is  decreased  by  1  are  given  by 


Xt  -  «J*i  ; i-1,2, . . .  ,d. 


(6-31) 


As  it  was  stated  earlier,  assuming  the  noise  components  to  be  stationary 
and  statically  independent,  it  can  be  shown  that  v^n  has  a  similar 
form;  i.e, 
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(6-32) 


(  ElSfllfl<n-1>BHA1HA1B'r<1-1)S]  ;i#n 

1  ElSH'fB<n-1)BHA1HA1Bf<i“1)S]  +  (m-d)ff2  ;i-n 


;i*n 
(6-33) 
;  i-n 

where 

Spq*E(SqSp] 

apq-ajap 

P'a  .mzdeJ(i-l)(^p  -  V  . 

M  i-1 

Let  then  P  and  Q  be  the  matrices 

P  -  ?i  -  (m-d)ff2  I 

Q  ■  Qj_  -  (m-d)o2  Ix.  (6-34) 

The  matrix  pencil  P-AN  becomes 

P-AQ  .  (P1-(m-d)o2I)-5(Q1-(m-d)ff2I1). 

It  can  be  shovn  that  P  and  Q  cam  be  decomposed  as 

p  -  u'v'u'8 

Q  -  U'v'l^u'8,  (6-35) 


Z  Z  SpqapqF'  e-j(n-l>*q 
q-i  p-i 

Z  Z  SpqapqF'pqe-3(n-l)+q  ej(i-D+p  +  (m-d)*2 
,o-l  P-1 
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vhere  f  was  defined  previously  and  u'  and  v'  are  the  matrices 


1  1  ...  1 

eJ*!  eJ^2  ...  e^d 


U' 


;j<d-i)*!  ;j(d-D*2 


^l(d-l)*d 


sllallF'll  s12a12F'l2 
s21a21F  21  s22a22F  22 


v' 


sdladlF  dl  sd2ad2F  d2 


sldaldF'ld 
s2da2dF  2d 


sddaddFdd 


The  matrix  pencil  P-SQ  then  becomes 

P-8Q  -(u'v'u'^-SOjVlflu'8) 

-U'v'd-STHju'8,  (6-36) 

which  satisfies  the  requirements  of  the  pencil  theorem  provided  all 
are  distincts.  Hence,  the  rank  reducing  numbers  of  the  pencil  P-5Q  are 
given  by 

i*  .  eJ*i  -  eJTwi  ;  i-1,2,.  .  .,d.  (6-37) 

Equation  (6-37)  together  with  equation  (6-31)  allows  us  to  estimate  si¬ 
multaneously  the  angular  frequencies  and  the  angles  of  arrival  of  the 
sources;  i.e, 
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■  -jln(Sj)/T 

©i  -  sin-1{(-jln(Xi))/(«1  D)} ;  i-1,2,.  •  -,d.  (6-38) 

6.3  COMPUTER  SIMULATION 

The  model  used  in  the  simulation  consisted  of  two  coherent 
sources  (d»2)  incident  on  a  linear  array  of  eight  uniformly  spaced 
sensors  (m-8).  For  convenience,  the  sensors  are  assumed  to  be  omnidirec¬ 
tional.  The  noise  was  simulated  as  white  Gaussian  with  zero-mean  and 
variance  o2  ■  1.  The  complex  envelopes  of  the  signals  emitted  by  the 
sources  are  assumed  to  be  constant  and  equal;  i.e, 
si(t)  •  s2(t)  *  s. 

The  sources  were  assumed  to  be  located  at  0^»180  and  02-22°  with  center 
frequencies  given  by  u^«0.2x2n  rd/s  and  »2«0.25x2n  rd/s  respectively.  D 
and  T  were  assumed  to  be  equal  to  6  108  meters  and  1  second  respective¬ 
ly. 

With  the  above  definition,  the  signal  to  noise  ratio  is  given  by 

2  Is  I2  |2 

SNR  .  -  -  2  s  2 

o2 

The  cases  considered  in  this  simulation  are  shown  in  table  (6-1) 
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Table  (6-1) 


SNR  | 

M 

30  dB  | 

22.36 

10  dB  | 

2.24 

Equation  (6-16)  vas  not  used  because  it  assumes  that  both  magnitudes  of 
X  and  6  are  unity.  Because  this  vas  not  the  case  in  actual  practice  due 
to  numerical  inaccuracies,  the  method  did  not  perform  veil.  One  vay  to 
solve  the  problem  is  to  use.  the  technique  as  in  chapter  4. 

Let  X  -  a  +jb  and  6  ■  c  +jd  vhere  a,  b,  c  and  d  are  real  numbers.  Using 
polar  form,  ve  have 

X  ■  {a^+b^)^exp{j tan“*(b/a)} 

and 

S  -  (c^+d2)^exp{j tan-l(d/c)} 


Prom  equations  (6-12)  and  (6-15),  X»exp{-j$)  and  S»exp{-j^}.  Thus,  ig¬ 
noring  the  magnitudes  of  X  and  & 

♦  -  -tan~l(b/a) 
and 


♦  »  -tan“l(d/c) . 


The  results  of  the  simulation  are  shovn  in  tables  (6-2)  and  (6-3). 


Table  (6-2) 


SNR 

1  “1 

«2  1 

©1 

©2 

30  dB 

|  0. 2001x2 n 

0. 2495x6n  | 

18.0012 

|  22.0233 

10  dB 

|  0. 1925x2 n 

0. 2448x6  It  | 

18.5130 

|  22.7785 

Mean  of  0^  and 

(  500  snapshots/run  ,  10  runs) 


Table  (6-3) 


SNR  |  Wj  u>2 

1  ©1 

©2 

30  dB  |  0.3687  10‘4  |  0.0546  10"4 

|  0.0111  I 

0.0045 

10  dB  |  0.0033  |  0.0004 

1  1.0504  I 

0.4662 

Variance  of  and 
(  500  snapshots/run  ,  10  runs) 

Observe  that  the  bias  in  the  estimates  is  very  small  at  both  high  and 
lov  signal  to  noise  ratio.  The  technique  can  thus,  be  classified  as 
giving  very  good  results. 
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CHAPTER  7 

I 

CONCLUSION 

In  this  study  the  objective  vas  to  extend  the  matrix  pencil  ap- 
'  proach  proposed  by  H.  OUIBRAHIM  (1].  The  flexibility,  the  eff ictiveness 

and  the  ease  of  use  of  the  method  are  once  again  shovn.  A  generalization 
the  method  to  arbitrary  but  identical  beam  patterns  was  derived  in 
chapter  3.  This  was  done  to  prove  that  the  choice  of  omnidirectional 
sensors  is  not  the  only  choice  possible.  In  chapter  4  we  used  different 
windows  to  generate  the  sequences  needed  in  the  formulation  of  the 
matrix  pencil.  It  was  shown  that  the  rectangular  window  performed  best 
compared  to  the  Hanning,  Hamming  and  Blackman  windows.  Errors  arise  when 
the  assumption  of  uniformly  spaced  elements  is  not  adhered  to  in  prac- 

I 

tice.  For  this  reason,  a  perturbation  analysis  due  to  sensor  spacing  was 
performed  in  chapter  5.  We  made  use  of  the  chordal  metric  introduced  by 
Stewart  113].  The  chordal  metric  provides  a  good  measure  between  the 
perturbed  and  the  unperturbed  eigenvalues.  It  was  shown  in  a  computer 
simulation  that  a  useful  bound  vas  found  for  the  chordal  metric.  In 
chapter  6  a  new  technique  was  introduced  to  simultaneously  estimate  the 
angular  frequencies  of  the  signals  emitted  by  the  sources  and  the  direc¬ 
tions  of  arrival  of  these  sources. 
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